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Abstract 

The  scattering  matrix  of  a 36O  degree  circumferential  slot 
in  an  infinitely  long  circular  waveguide  with  an  incident  dominant 
circular-electric  wave  is  obtained  by  a variational  principle* 

The  theory  which  should  be  very  good  for  small  gaps  is  shown  to  be 
in  good  agreement  with  the  obtained  experimental  results* 


I.  Introduction 

Consider  a circumferential  gap  in  an  infinitely  long,  per- 
fectly conducting  circular  waveguide  as  shown  in  Fig*  !<.  The 
gap  can  be  considered  to  be  a region  of  free-space  coupling  two 
serai-infinite  circular  waveguides.  Assuming  that  the  gap  width 
is  small  compared  to  the  guide  wavelength,  the  basic  problem  is 
to  determine  the  effect  of  the  gap  on  an  incident  TE01  wave* 

The  energy  in  the  IEq^  wave  incident  on  the  gap  is  partially 
reflected  and  partially  transmitted.  Some  of  the  energy  is  cou- 
pled through  the  gap  to  the  free  space  surrounding  the  waveguide* 
Since  the  power  relationships  are  the  essential  quantities  of 
interest,  the  problem  is  formulated  in  terms  of  the  scattering 
matrix  of  a waveguide  junction. 


•-  * ’.ki.-xm 
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Intuitively,  the  small  gap  should  have  little  effect  on 
the  propagation  of  an  Incident  circular-electric  wave.  The 
circular-electric  modes  have  the  unique  property  that  the  wall 
surface  currents  are  purely  transverse.  These  transverse  cur- 
rents do  not  tend  to  charge  the  metallic  boundary  of  the  circum- 
ferential gap$  hence  little  energy  is  expected  to  radiate  into 
space.  On  the  other  hand,  all  TM  modes  are  characterized  by 
solely  longitudinal  currents  on  the  guide  wall,  and  all  TE  modes 
other  than  the  circular-electric  type  have  longitudinal  components 
of  current.  It  can  be  concluded  that  the  modal  energy  of  all  ex- 
traneous modes  will  be  coupled  strongly  to  the  surrounding  medium. 

The  problem  is  formulated  in  terms  of  an  integral  equation 
for  the  axial  component  of  the  magnetic  field  on  the  gap  surface. 

A variational  principle  is  applied  and  a parameter  A is  introduced 
representing  the  stationary  parameter  in  the  variational  formulation 
of  the  integral  equation.  For  simplification,  even  and  odd  excita- 
tions are  used  and  the  elements  of  the  scattering  matrix  are  ex- 
pressed In  terms  of  the  appropriate  A*s.  The  restriction  imposed 
by  the  theory  is  that  the  gap  width  must  be  small  compared  with 
wavelength. 


2.  Formulation  of  the  Problem. 

Since  the  following  analysis  pertains  to  circularly  symmetrical 
waveguide  modes,  the  electromagnetic  field  components  can  be  de- 
rived from  a scalar  function  which  is  independent  of  the  fb-  coor- 
dinate. Assuming  a time  dependence  exp  -iwt,  the  circular-electric 
modes  are  obtained  by  solving  the  equation 


L B^(r,z)  * 0 


(1) 


where  the  operator 
given  by 

L 


L acting  on  the  ^-components  of  the  field  is 


I 


l 


FIG.  i 


A CIRCUMFERENTIAL  GAP  IN  AN  INFINITELY 
LONG  CIRCULAR  WAVEGUIDE 


TR166 


-3- 


The  other  field  components  can  be  expressed  in  terms  of 
E^  by  5 

Hr(r*z)  * ir  h Vr>z) 

Hz(r’z)  ' j|i  r ir  rVr,z)  • (3) 

E„  * E„  “ Hv  * 0 
T Z P 

The  walls  of  the  cylinder  are  assumed  to  be  perfectly  con- 
ducting! therefore,  the  boundary  conditions  at  the  wall  require 
that  the  tangential  components  of  the  electric  field  vanish. 
Hence,  the  boundary  condition  is 

E^(r,z)l  = 0 (4) 

J r*a 

A circumferential  gap,  of  width  2g,  centered  at  z c 0 is 
introduced  in  the  infinite  guide  as  indicated  in  Fig.  2.  The 
region  for  r < a *ill  be  designated  by  superscript  I and  the 
region  exterior  to  the  guide,  r > a,  by  superscript  IT. 

In  region  II  (r  > a)  the  field  can  be  expressed  in  terms  of 
a Green* s function  and  the  field  distribution  on  the  surface  of 
the  gap.  The  Green's  function  is  defined  by  the  equation 


L 011  (r,r'  ,z-z' ) - - »(r-r‘Wz-tU  r>a 


G11  (r,r’,z-z 


■L 


= 0 


(5) 


By  making  use  of  Green's  second  theorem,  an  integral  represen- 
tation of  the  external  field  can  be  found  in  the  form 


(r*z)  = J°  [h  5FT  r* 


G11  (r,r',z-z') 


? (z')dz' 

r‘=a  (6) 


where 


% (z‘)  c [s^(r  ,z')] 


r=a 


i 

i 

i 

! 
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While  an  explicit  form  for  this  Green’s  function  can  be 
obtained  if  desired,  it  is  not  necessary.,  Actually,  it  is 
simpler  to  determine  the  function  occurring  in  the  integral  of 
(6).  It  is  apparent  from  (6)  that  the  integrand  must  satisfy 
the  two  conditions 


L 'h  sfr  r'  GTI(r,r',z-z»)  | 

L “*  r ! *a 


[* 


— r ' GIT(r 


dr  ’ 


,r • ,z-z ' ) 1 

J «t 


0 r > a 


5(z~z ' ) 


(7) 


(8) 


r "*a 
r “a 


as  well  as  the  Sommerfeld  radiation  condition.,  By  the  use  of 
the  Fourier  transform  one  obtains  an  expression  of  the  form1 

"l  3 i II  1 i T®  it(z-z’)  H^1)  (lAc^-  5 2 r) 

["rT  5rT  r (r’r,*M,)J  8 jn  n/>j  gya)dg- 


where 


Vk2  - 5 2 « £V^  2 - k2  for 


(9) 


It  is  immediately  apparent  that  such  a choice  satisfies  (7) 
and  (8)  as  well  as  the  radiation  conditions „ Since  the  Green's 
function  is  unique,  (9)  is,  therefore,  the  desired  function., 

A similar  analysis  can  also  be  carried  out  for  the  region  I 
(r  > a)0  Since  there  are  fields  propagating  inside  the  pipe, 
the  resulting  expressions  for  the  fields  include  an  additional 
term  representing  the  incident  wave*  In  this  case  one  obtains 


B^(r,z)  = 8 


I 

<6  Inc 


(r,z) 


Performing  an  analysis 
possible  to  show  that  ‘ 


*(r,r'  ,z-z'  )J 


g ( z ' )dz’ 
r ' =a 


(10) 

similar  to  that  for  region  II,  it  is 


REGION  H 

EXTERIOR  TO  GUIDE 
AND  GAP  SURFACE 


INFINITELY  LONG 
CIRCULAR  WAVEGUIDE 


Z = +g 


Z = 0 


REGION  X 

INTERIOR  TO  GUIDE 
AND  GAP  SURFACE 


FIG.  2 ORIENTATION  OF  CYLINDRICAL  COORDINATE  SYSTEM 


r_k.  »_ 

L rs  <Sr 


G*(r,r‘  ,z.-z’ ) 


] 


r ,=a 
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iS(z-z»)  J1(Vk5-?2  p) 

2 JlCTO^ 

(11) 

Assuming  only  a single  mode  propagates  in  the  pipe,  the 
Bessel  function  in  the  denominator  of  (11)  vanishes  at  two  places 
(+  Y][)  on  the  real  axis,  and  at  an  infinite  set  of  points  (+  iyn) 
along  the  imaginary  axis.  The  contour,  C,  used  in  (11)  is  thau 
indicated  in  Pig.  (3).  The  contour  is  chosen  to  pass  around  the 
singularities  in  the  manner  shown  so  that  the  Green’ is  function 
represents  only  waves  propagating  away  from  the  junction,. 

Thus,  the  entire  electromagnetic  field  has  been  expressed  in 
terms  of  the  field  on  the  surface  of  the  gap. 
determine  integral  equations  for  these  fields, 
by  matching  the  magnetic  fields  across  the  gap. 


It  remains  now  to 
This  can  be  done 


d5 


In  particular,  from  (9)  the  magnetic  field  in  the  gap  surface 
can  be  shown  (except  for  dimensional  constants)  to  be 


where 


[il 

l_r  d r 


rE^I(r,z) 


r-a 


i 


KII(z-z')  £(z')dz'  (12) 


KTI(z-z ' ) 


= JL . 

2nr 


oo 

/ 

W-00 


i^(z-z') 
e H(5) 


(13) 


and 


H(S) 


H a) 


E. 


(1) 


(Vk2-S‘ 


a) 


(14) 


Similarly,  from  (10) 
Hz(a.z>  *[r  & T EJ<r*z>] 


" [ r Jr  r E^Inc.^r’z^l 

r=a  r=a 


t 
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+ J*  KT(z-z» ) \ ( z ’ ) d z ’ 


KT(2-2.)  « £ f 


l$(z-z * ) 


e 


J(5)  = 


_ o 


C 

0/kJ-?2 


a) 


J^Vk2-^2  a) 


(15) 


(16) 


(17) 


Hz  <aiz)  = ^(a*2)  -g  < z<g  , 
and,  therefore,  one  obtains  the  integral  equation 


r & r EJlnc.(r.z> 

‘ f 

KIT(z-z»)  - KT(z-z») 

r= a «-g 

£ (z')dz'  -g<  z < g 

(18) 


3 • Variational  Formulation. 

It  Is  now  convenient  to  separate  KX(z-z')  Into  two  terms 
follows : 


as 


Kr(z-Z')  . ftz-z’)  f t 1?(Z"Z,) 


(19) 


where 


rl/  % 1 i5(z-z») 

^(z-z')=-^  / e J($)d$ 


(20) 


■00 


i.e. , integrals  along  the  real  axis  and  similar  in  form  to  those 
of  (13).  The  remaining  integrals  along  the  contour,  C,  indicated 
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in  Pigo  4 can  be  evaluated  by  means  of  residues,. 
The  integral  equation  then  assumes  the  form 


iY,a 


■ ao  r 

2 r fr  r EJlnc.(r’z^ 
1 L 


] +/ 


cosy^z-z*  )£(z  ‘ )dz ' 


K(z-z  * )|(z 5 )dz ' 


where 


Y a-5 

K(z-z')  = -Jy-  KIX(z-z')  - Kr(z-z') 

ipj  L J 

and,  the  constants  P19  Y1  are  defined  by 

* ri  *Vk2  - <^)2 

The  fields  incident  in  the  waveguide  can  be  expressed  in 
the  form 


!ilnc.“Nz)  ■ 


_fl_  Jl(plf>  , . 
, -2  "7700  a(z) 


iY^  V*V 

Their  z-dependence^  a(z),  is,  in  general,  of  the  form 

, e iYlZ  -iY  Z 

a(z)  = Ae  x + Be  1 

i.e.,  fields  incident  on  the  gap  from  either  the  right  or  left, 
It  is  convenient  to  specialize  this  into  even  or  odd  incident 


a^Cz)  = cos  Y^z 


OqCz)  = sin 


even 


For  the  even  case,  £e(z)  = {-  (-z)  and  hence. 


^ V'“  * 


;s'-^(z-z 1 )£eCz  ')dz'  = 


ac(«) 


/ °e<2,;' 


% (z* )dz* 


Similarly,  for  the  case  of  odd  excitation  £ (z)  = - £Q(-z) 
and  hence. 


cosY. 


re 

(z-z* )|0(z' )dz'  = aQ(z)  / a0(z* )^Q(z ' )dz ' 


Therefore,  for  either  even  or  odd  excitation,  1.23)  can  be  inserted 
in  the  integral  equation  (21)  which  then  becomes  s 


a(z)  1 + / a(z' )|(z' )dz*  = / K(z-z’ )|(z' )dz • 


• (24) 


Multiplying  (24)  by  £ (z),  integrating  from  -g  to  g and 
dividing  by 


a(z)|(z)dz  yields  the  variational 


principle 


%* 


*/:■ 


(z)|(z)dz 


rj. 


a(  z)j*  (z)dz 


K(z-z,)j'(z)j'(zJ)dzdz' 


a(z)^(z)d; 


(25) 


(25)  can  be  shown  to  have  all  the  usual  variational  properties, 
i.e.,  A is  stationary  with  respect  to  small  changes  of|(z)  about 
its  correct  value  and  A is  independent  of  the  overall  amplitude 


"T*r* 
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of  £<z) 


Thus,  it  is  possible  to  insert,  an  approximate  value  for 
£ (2)  in  (25)  and  obtain  relatively  more  accurate  values  of  A, 


As  will  be  seen  later,  the  various  A‘s  can  be  related  to  the 


scattering  matrix  of  the  Junction* 


Referring  to  equations  (13),  (14),  (17),  (20),  and  (24),  it 
is  seen  that  K(z -z5)  can  be  expressed  in  the  form 


'>00  iS(z-z') 

K(z-z')  = I e *($)d£ 


(26) 


/-CD 


where,  using  the  Wronskian  for  Bessel  Functions, 

2 


*(5)  = a)H^W  a) 


(27) 


The  variational  principles  can  then  be  written  in  the  form 


L = I 
'-00 


(28) 


where 


cos  $z  |(z)dz 


1.(6)  ■ 


e 


(29) 


cosy-ji  z £ (z)dz 


'-g 


sin  £ a | (z)dz 


§0<*;>  - ^ 


r sinY^^z  ^z)dz 

J-g 


It  remains  now  to  choose  trial  functions  forl-(z)  and  insert 
them  into  (28),, 


n 


•10- 
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Irnsider  the  follorlng  choices 

|eCz)  = Vg2  -z2  !0(z)  = z Yp~~^P  VIC 

It  Is  apparent  that  these  functions  satisfy  the  symmetry  conditions 
and  have  the  correct  static  dependence  on  the  °:ge  the  gao. 

small  gaos  they  should,  therefore,  be  good  approximations  to 
the  correct  fields.  Inserting  these  into  (29)  yields 


*e<«  = 


sji<ng) 


2 


(fcg) 

(^g)2 


5 - 


VY1K) 

(y^e' 

(3D 


Inserting  (31)  into  (28),  it  can  be  shown  (see  Appendix)  that 


= 


Yla 


( K.9.).  + (~)2j 

4 TT  Vg'  1 


£ ^ — 
o Y 


laPi 


(&)4 

16  n vg' 


It  is  to  be  noted  that  these  results  are  only  valid  for  small 
2 2 

gaps,  i.e.,  g /\  <<  1.  Results  for  larger  gaps  would  be  difficult 

to  obtain  theoretically  because  of  the  lack  of  satisfactory  trial 
functions  for  £(z).  Moveover,  it  is  virtually  impossible  to  make 
satisfactory  experimental  measurements  for  large  gaps  since  the 
interaction  of  the  gap  with  the  "free  space"  region  outside  the 
guide  is  of  prime  Importance  and  a "free  space"  region  is  difficult 
to  duplicate  experimentally. 


I 
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4.  Transmission-Line  Representation 


It  is  now  desirable  to  express  the  scattering  matrix  ele- 
ments in  terms  of  and  AQ.  This  can  be  accomplished  by  ex- 
amining the  fields  inside  the  waveguide  at  a terminal  plane  lo- 
cated a great  distance  from  the  gap  discontinuity.  The  evalu- 
ation of  the  scattering  (or  impedance)  matrix  requires  the  knowl- 
edge of  the  explicit  asymptotic  expression  of  the  Green’s  function 
of  equation  (11).  This  expression  is  obtained  in  the  appendix. 
Substitution  of  the  asymptotic  form  into  equation  (10)  yields 


EI(r  z)  = E1  (r  z)  + J1(P}  e '1Y1Z  f 

E^(r,z;  E^inc.(r^z;  Yja2  e J 

& 


iY-j  z ’ 

e 1 £(z')dz’ 


(33) 


for  the  electric  field  in  the  far  zone. 

It  then  follows  that  the  z dependence  is  of  the  form 

— 1 y z i Y z ’ 

E^(r,z)^a(z)  + e 1 / e 1 ^(z’Jdz*  (34) 

J-  g 

Now  for  the  even  case  (a(z)  = cosY^z)  this  expression  becomes 


T iY^  z 

E^(r,z)  ~ e + 


J-g 


1 + 2 cosY^z ’ |(z')dz 


e 


-iY^z 


(35) 


and  the  even  reflection  coefficient  P^  can  be  defined  as 

e 


P0  = 1 + 2 

e • 


& 


f-g 


cosy^z’  ^(z’)dz' 


(36) 


Similarly,  for  odd  excitation  a (z ) = sin  y^z,  the  odd  reflection 
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coefficient  ro  can  be  expressed  as 


O 

no  = - 1 - 2 siny^z*  £(z’)dz' 


(37) 


The  integrals  of  equations  (36)  and  (37)  can  be  expressed  in 
terms  of  A from  equation  (25),  resulting  in 


A.  - i 

P « -T-T 
1 e L T-  1 

O 

A-  - 1 


(38) 


f’o-tJS 


The  elements  of  the  scattering  matrix  are  then  given  by 


S11  “ S22 


P + n 

-A,..  1 9.. 


1 + A„  1 + A 
0 e 


*12 


r«  ~ r0  . x . 1 1 

2 ttfa;  tts; 


(39) 


5o  Experimental  and  Theoretical  Results 

Measurements  of  the  characteristics  of  the  circumferential 
gap  were  performed  in  a 5» 755-in.  i.d.  brass  tubing  at  a free- 
space  wavelength  of  10.000  cm.  Deschamps'  method  was  used  to 
determine  the  elements  of  the  scattering  matrix  of  the  circum- 
ferential gap.  In  comparing  the  theoretical  and  experimental 
results  it  was  convenient  to  examine  separately  the  absolute 
values  and  arguments  of  the  scattering  matrix  elements.  This 
method  of  presentation  of  data  was  chosen  because  the  experi- 
mental values  of  amplitude  and  phase  were  determined  separately 

p 

by  Deschamps’  method. 


AMPLITUDE  OF  S„  | | H PHASE  OF 


FIG.  5 MEASURED  AND  CALCULATED  VALUES  OF  S 


•» 


> y**  mv i cmmmrnmu 
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The  absolute  values  of  S 22 ? an<*  are  ^en  §*ven  by 


lsul 


= Is 


22 


Aa  - A 
e o 


1 + A„  1 + A 
o e 


(40) 


’12 


Ka 

- 1 

1 * Ao 

1 

+ AJ 

and  their  arguments  are  found  to  be 


Arg  = Arg  S22  = Arg(A0  - Afi)  - Arg(l  + AQ)  - Arg(l  + A0) 

Arg  S12  = Arg(A0Ao  - 1)  - Arg(l  + AQ)  - Arg(l  + A0)  (4-1) 

The  experimental  values  are  plotted  on  the  same  figure  as  the 
theoretically  determined  curves  in  Figs.  5 and  6.*  It  is  observed 
that  the  experimental  value  of  (S'^j  and  |S12|  are  in  excellent 
agreement  frith  the  theory  for  values  of  2g  up  to  1.000  in.  corre- 
sponding to  £ * 0.127.  It  is  recalled  that  the  theory  is  only 

A.  p 

claimed  to  be  valid  for  values  of  «1.  As  seen  from  the 

phase  characteristics  of  S^,  the  experimental  values  lie  slightly 
above  the  theoretical  curve  for  small  values  of  2g.  It  is  observed 
in  Fig.  6 that  there  is  again  excellent  agreement  between  experiment 
and  theory  for  values  of  Arg  S^2  such  that  2g  1.000. 


t 


I • 


A detailed  description  of  the  experimental  techniques  as  well 
as  a more  complete  discussion  of  the  experimental  results  is 
given  in  Cruft  Laboratory  Technical  Report  No.  167.  It  should 
also  be  noted  that  in  calculating  the  theoretical  values  of  the 
scattering  matrix  elements  from  equation  (32)  is  substituted 
for  + i to  conform  to  the  convention  established  in  measurement# 
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Appendix  A 

Evaluation  of  the  Variational  Parameters  ^e  Mi  4^ 


The  determination  of  the  elements  of  the  scattering  matrix 
of  a circumferential  gap  requires  the  evaluation  of 


A 


o 


(A-l) 


kince  the  assumed  electric  field  distribution  in  the  gap  is 
a good  approximation  only  for  small  values  of  Y^g,  the  parameters 
5 (5)  and  * (§)  can  be  expressed  in  the  approximate  form 

© O 


*e(S>  = 


__  V^Sg)  c 2J1(^g) 


S^O^g)  £g 


r1g«  l 


(A-2) 


• 0C5)  = 


2 (t  . Jo^g) 

2 ^ x(^g ) ” 


i _a_ 


2jlCYXg)  J0(Y1g)  Yxg 
(Y^)2  Y1B 

and  ir(&)  is  given  by 

n$)  = 


C£g) 


fcg 

J 

ng  « 1 


V‘ 


a) 


(A-3) 


It  is  convenient  to  determine  separately  the  real  and  imagi- 
nary parts  of  Ag  and  AQ.  Proceeding  to  evaluate  R.P.  Ae  one 
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observes  that  the  integrand  is  real  only  for  £ <k.  Therefore, 
assuming  kg  <<1  as  well  as  Y^g  <<19  so  that  J^(^g)  = ^ for 
0 < t,  < k,  it  follows  that 


2V  a2 


R.P.  A = — ?r 
e __<■ 


^ P1  4>  ^ a)  + Y^( 


(A-4) 


V - s a) 


2.  1 


2V  i r dt 

n^Pl  1 |H^(/TT7ka)|2 


For  t<<i,'  the  denominator  is  large,  and  hence  there  is 
little  contribution.  On  the  other  hand  for  t near  zero,  and  with 
ka  equal  to  4.592,  the  Hankel  function  is  already  well  approximated 
by  its  asymptotic  form 


Vl-22  (kal  * 


=r  1 + 

C?  L 


8(l-t2)(ka)‘ 


(A-5) 


Inserting  (A-5)  Into  (A-4)  results  in 


R.P.  A = 
e 


i Yi*^a3  f 

Jo  1 + r 


(A-6) 


y ■ 

*> 
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where  the  quantity  6 is  defined  as 


8(ka)‘ 

Equation  (A-6)  can  be  written  as 


■ . vv 

R.P.  A„  = — 

6 4P  2 
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n J 5 + (1-t2) 
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Now  since  5 <<  1,  it  can  be  neglected  in  evaluating  the  remaining 
integral  of  (A-7).  One  then  obtains 


2 R 

■ ••nil  r-  -»  Y.k  a 

R.P.  A * — — 5-  ll  - 26  I = — — 5- 


Y,k2a3 


4p; 


4P 


(A-8) 


1 "rl 

To  evaluate  the  imaginary  part  of  A , examine  the  integral 


Y.a2  ff° 

A c T D i / 
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(A-9) 


Noting  that  the  integrand  is  imaginary  for  values  of  § > k,  and  that 
the  major  contribution  to  the  integral  occurs  for  £ >>  k,  one  can 
write 


jT6[£ 


(A-10) 


where  the  asymptotic  expression 

Jx(Vk2-^2  a )E[1)0fk2-^2  a)  * JjU&OH^CiSa) 


(A-ll) 


niga 


Ua|  > > 1 


has  been  used.  It  then  follows  that* 


I.P.  A_  = 
e 


8r,a  . , f 

4 — | 
>*1  * Jo 


[J1c^g)j 


d$  = 


♦Watson,  A Treatise  on  the  Theory  of  Bessel  Functions.  Cambridge 
University  Press,  London,  i94?,  p.  397. 
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In  a similar  fashion  the  real  part  of  A is  given 


by 


R.p.  ao  = 


*1»2,  8 2 -® 


R-p-  f 

n pi  Ylg  L 


“till2  Ji(58)  - s 


2 


g 


00  a )Hj1)(Vk2-^2  a) 


d£ 


(A-13) 

It  is  seen  again  that  the  integrand  has  a real  part  when  k < $ . 

gain  assuming  kg  <<  1,  and  using  the  leading  term  in  the  extension 
of  the  numerator,  one  obtains 


R.P.  A. 


~~2T 

"PlY 


52 


i UO 


H^(Vk2  - 


r a) 


d5 


Again  making  the  substitution  t = 5/k, 


n n 2(ka)3 

R.P.  A_  = —*  

0 ,2  2. 


t2dt 


n P 1 V J0  |H^(l/l-t*  ka)  | 2 

Introducing  the  asymptotic  expression  I (Vl-t2  ka)| 
this  becomes  * 1 I 


(A-l4  ) 
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Finally,  by  using  an  argument  similar  to  that  which  led  from 

(A -9)  to  (A-10),  the  imaginary  part  of  A is  seen  to  be  of 
the  form 


where 


x-p-  A0  » -P-cf)4  I 

nPlYla  * 


I \&  ^ ^ 


By  the  substitution  t * £g. 


J Ji(t)  - $ ♦ ^r- 


This  integral  can  be  expressed  as 


f j?<t> 
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Therefore 


I.P.  a = -22-  (ft) 


”pjr1a  * 


Finally,  the  resultant  expressions  for  the  variational 
parameters  are  given  by 


i 7|  + 1 n (t)2J  • 

+ 1 
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Appendix  B 


Evaluation  of  the  Explicit 


Expression  of 


the  Green's  Function  of  Equation  (11) 

It  is  desired  to  derive  an  explicit  form  for  the  Green's 
function  of  equation  (11)  in  order  to  evaluate  the  far-zone 
electric  field  in  region  I.  The  expression  for  the  Green's 
function  is  given  by  /~5 — ~ 


* rt  Sp*  r,Gl(r9r'  <z~z' ) 


= 


+i£(z-z')  j m m 


J,(yk  -5'a) 


(B-l) 


where  the  contour  C is  chosen  as  in  Fige  3 


It  is  noted  that  the  integrand  is  analytic  except  for  simple 
poles  at  £ « + Y^  and  £ ■ + IYn  (n>  1).  Since  the  result  is 
required  lor  z «z'  it  is'*  appropriate  to  close  the  contour  with 
a large  semi-circle  in  the  lower  half- plane . It  is  apparent 
that  the  integral  along  the  semi-circle  vanishes  as  its  radius 
tends  to  infinity*  Therefore,  by  the  theory  of  residues,  and 
noting  that  only  the  poles  % = -Y^  and  % **  - IY  (n  > 1)  are 
enclosed,  one  obtains 

" JT  syr  r'GI(r,r',z-z‘)  = §“■  Residue  at 

L Jr'=a  *Y1  (B-2) 

$*  “ iYn  (n^l) 

j 2_z » ) 

•Since  e - J represents  a rapidly  attenuated  function  when 

z «z'  and  £ is  negative  imaginary,  all  higher-order  terms  are 
negligible  compared  to  the  dominant-mode  term. 

Evaluating  the  residue  corresponding  to£=  "Y^,  une  obtains 
the  desired  expression 


iY  (n  > 1)  are 
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